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Abstract 

We present recent evolutions of the detailed opacity code SCO-RCG which combines statistical mod¬ 
elings of levels and lines with fine-structure calculations. The code now includes the Partially Resolved 
Transition Array model, which allows one to replace a complex transition array by a small-scale detailed 
calculation preserving energy and variance of the genuine transition array and yielding improved high- 
order moments. An approximate method for studying the impact of strong magnetic field on opacity 
and emissivity was also recently implemented. The Zeeman line profile is modeled by fourth-order Gram- 
Charlier expansion series, which is a Gaussian multiplied by a linear combination of Hermite polynomials. 
Electron collisional line broadening is often modeled by a Lorentzian function and one has to calculate 
the convolution of a Lorentzian with Gram-Charlier distribution for a huge number of spectral lines. 
Since the numerical cost of the direct convolution would be prohibitive, we propose, in order to obtain 
the resulting profile, a fast and precise algorithm, relying on a representation of the Gaussian by cubic 
splines. 
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1 Introduction 


When atomic spectral lines coalesce into broad unresolved patterns due to physical broadening mecha¬ 
nisms (Stark effect, auto-ionization, etc.), they can be handled by so-called statistical methods (Bauche et 
al., 1988). Global characteristics - average energy, variance, asymmetry and sharpness - of level-energy, 
absorption or emission spectra can be useful for their analysis and the investigation of their regulari¬ 
ties (Bauche & Bauche-Arnoult, 1987 & 1990). Systematic studies of these average characteristics for 
transition arrays and applications to the interpretation of experimental spectra of high-temperature 
plasmas were initiated in (Moszkowski 1962, Bauche et al., 1979). The elaboration of the general group- 
diagrammatic summation method (Ginocchio, 1973) and its realization in computer codes (Kucas & 
Karazija, 1993 & 1995, Kucas et al., 2005, Karazija & Kucas, 2013) opened up new possibilities for 
the use of global properties in atomic spectroscopy. On the other hand, some transition arrays exhibit 
a small number of lines that must be taken into account individually. Those lines are important for 
plasma diagnostics, interpretation of spectroscopy experiments and for calculating the Rosseland mean 
kr, important for radiation transport, and defined as 
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hv being the incident photon energy, n{hi') the opacity including stimulated emission, T the temperature 
and Rt(^J') Planck’s distribution function. The Rosseland mean is very sensitive to the gaps between 
lines in the spectrum. 

These are the reasons why we developed the hybrid opacity code SCO-RCG (Porcherot et al., 2011), 
which combines statistical methods and fine-structure calculations, assuming local thermodynamic equi¬ 
librium. The main features of the code are described in section [2j the extension to the hybrid approach 
of the Partially Resolved Transition Array (PRTA) model (Iglesias & Sonnad, 2012), which enables one 
to replace many statistical transition arrays by small-scale DLA (Detailed Line Accounting) calculations, 
is presented in section [3] and comparisons with experimental spectra are shown and discussed in section 
U In section [5j an approximate modeling of Zeeman effect is proposed together with a fast numerical 
algorithm for the convolution of a Lorentzian function with Gram-Charlier expansion series, based on a 
cubic-spline representation of the Gaussian. 


2 Description of the code and effect of detailed lines 

In order to decide, for each transition array, whether a detailed treatment of lines is necessary or not and 
to determine the validity of statistical methods, the SCO-RCG code uses criteria to quantify the porosity 
(localized absence of lines) of transition arrays. The main quantity involved in the decision process is 
the ratio between the individual line width and the average energy gap between two lines in a transition 
array. Data required for the calculation of lines (Slater, spin-orbit and dipolar integrals) are provided by 
SCO (Superconfiguration Code for Opacity) code (Blenski et al., 2000), which takes into account plasma 
screening and density effects on the wave-functions. Then, level energies and lines are calculated by an 
adapted routine (RCG) of Cowan’s atomic-structure code (Cowan, 1981) performing the diagonalization 
of the Hamiltonian matrix. Transition arrays for which a DLA treatment is not required or impossible 
are described statistically, by UTA (Unresolved Transition Array, Bauche-Arnoult et al., 1979), SOSA 
(Spin-Orbit Split Array, Bauche-Arnoult et al., 1985) or STA (Super Transition Array, Bar-Shalom et 
al., 1989) formalisms used in SCO. In SCO-RCG, the orbitals are treated individually up to a certain 
limit, consistent with Inglis-Teller limit (Inglis & Teller, 1939), beyond which they are gathered in a 
single super-shell. The grouped orbitals are chosen so that they weakly interact with inner orbitals (this 
is why we sometimes name that super-shell “Rydberg” super-shell). The total opacity is the sum of 
photo-ionization, inverse Bremsstrahlung and Thomson scattering spectra calculated by SCO code and 
a photo-excitation spectrum in the form 
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where h is Planck’s constant, M the Avogadro number, eo the vacuum polarizability, m the electron 
mass, A the atomic number and c the speed of light. V is a probability, / an oscillator strength, ^(hv) 
a profile and the sum X —> X 1 runs over lines, UTA, SOSA or STAs of all ion charge states present in 
the plasma. Special care is taken to calculate appropriately the probability of X (which can be either a 
level aJ , a configuration C or a superconfiguration S ) because it can be the starting point for different 
transitions (DLA, UTA, SOSA, STA). In order to ensure the normalization of probabilities, we introduce 
three disjoint ensembles: V (detailed levels aJ), C (configurations C too complex to be detailed) and S 
(superconfigurations S that do not reduce to ordinary configurations). The total partition function then 
reads 


u t ot = u(v) + u(C) + u(S) with vncns = < 


(3) 




where H is the Hamiltonian, 


where each term is a trace over quantum states of the form Tr 

N the number operator, /x the chemical potential and /? = 1/ (fcsT). The probabilities of the different 
species of the ./V-electron ion are 


V a j = -l-(2J+l)e-K E ‘> J -» N \ (4) 

Utot 

for a level belonging to P, 

V c = jj- (5) 

t-Aot Tr . n 
olJSLC 

for a configuration that can be detailed, 

V C = jf—gc e~^ Ec ~^ ( 6 ) 

for a configuration that can not be detailed ( i.e . belonging to C) and 

V s = jy— ^2 9° e~^ Ec ~^ (7) 

Utot C6S 

for a superconfiguration. 

We can see in Fig. |T] that fine-structure calculations can have a strong impact on the Rosseland mean. 
The physical broadening mechanisms are the same for both calculations (statistical: SCO and detailed: 
SCO-RCG). The modelling of the (impact) collisional broadening relies on the Baranger formulation 
(Baranger, 1958) and expressions provided by Dimitrijevic and Konjevic (Dimitrijevic & Konjevic, 1987) 
corrected by inelastic Gaunt factors similar, for high energies, to the ones proposed by Griem (1962 & 
1968). Ionic Stark effect is treated in the quasi-static approximation following an approach proposed by 
Rozsnyai (1977), corrected in order to reproduce the exact second-order moment of the electric micro-field 
distribution in the framework of the OCP (One Component Plasma) model (Iglesias et ah, 1983). 

The code can be useful for astrophysical applications (Gilles et ah, 2011, Turck-Chieze et ah, 2011). 
Figures [21 [3] and U] represent the different contributions to opacity (DLA, statistical and PRTA) for an iron 
plasma in conditions corresponding to the boundary of the convective zone of the Sun, with a maximum 
imposed number -/V max of detailed lines per transition array equal to 100 and 800000 respectively. As 
expected, when N max increases, the statistical part becomes smaller. In Fig. [21 the detailed part is 
obviously not sufficiently larger than the statistical part around /u/=850 eV. 

As shown on figures 0 [5] and [3 the statistical calculation (SCO) may depart significantly from the 
detailed one (SCO-RCG), and the differences are essentially the signature of the porosity of transition 
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Figure 1: (Color online) Comparison between the SCO-RCG and full-statistical (SCO) calculations for 
an iron plasma at T=50 eV and p=10~ 3 g/cm 2 . The Rosseland mean is equal to 1691 cm 2 /g for the 
SCO calculation and to 1261 cm 2 /g for SCO-RCG. 

arrays. The density being quite low ( p= 4 10 -3 g.crn -3 ), the lines emerge clearly in the spectrum. These 
conditions are accessible to laser spectroscopy experiments (see Sec. G} and we can see that the opacity 
changes notably with temperature. Therefore, even if the quantity which is measured in absorption 
point-projection-spectroscopy experiments is not the opacity itself, but the transmission (see Sec. HI): one 
may expect to have a reliable idea of the plasma temperature during the measurement. 


3 Adaptation of the PRTA model to the hybrid approach 

In order to complement DLA efforts, the code was recently improved (Pain et ah, 2013a & 2015) with 
the PRTA (Partially Resolved Transition Array) model (Iglesias & Sonnad, 2012), which may replace the 
single feature of a UTA by a small-scale detailed transition array that conserves the known transition- 
array properties (energy and variance) and yields improved higher-order moments. In the PRTA approach, 
open subshells are split in two groups. The main group includes the active electrons and those electrons 
that couple strongly with the active ones. The other subshells are relegated to the secondary group. A 
small-scale DLA calculation is performed for the main group (assuming therefore that the subshells in 
the secondary group are closed) and a statistical approach for the secondary group assigns the missing 
UTA variance to the lines. In the case where the transition C —> C is a UTA that can be replaced by a 
PRTA (see Fig. 0, its contribution to the opacity is modified according to 

fC-yC' 'Pc 'P&J (8) 

a.J—tOi'J' 

where the sum runs over PRTA lines aJ —> a’J 1 between pseudo-levels of the reduced configurations, 
faj-ya’j' the corresponding oscillator strength and '&aj-yd'j> is the line profile augmented with the 
statistical width due to the other (non included) spectator subshells. The probability of the pseudo-level 
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Figure 2: (Color online) Different contributions to opacity calculated by SCO-RCG code for an iron 
plasma at T=193 eV and p=0.58 g.cm -3 (boundary of the convective zone of the Sun). The maximum 
number of lines potentially detailed per transition array is chosen equal to 100. 



Energy (eV) 


Figure 3: (Color online) Different contributions to opacity calculated by SCO-RCG code for an iron 
plasma at T=193 eV and p=0.58 g.cm -3 (boundary of the convective zone of the Sun). The maximum 
number of lines potentially detailed per transition array is chosen equal to 800000. 
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Figure 4: (Color online) Comparison between the full-statistical (SCO) spectrum and SCO-RCG around 
the maximum of the opacity bump in the conditions of Fig. [2] (boundary of the convective zone of the 
Sun). The maximum number of lines potentially detailed per transition array is chosen equal to 800000. 



Figure 5: (Color online) Iron opacity at T=15 eV and p= 4 10~ 3 g.cm -3 . Comparison between the 
full-statistical (SCO) and SCO-RCG calculations. The maximum number of lines potentially detailed 
per transition array is chosen equal to 800000. 
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Figure 6: (Color online) Iron opacity at T= 27 eV and p= 4 10 -3 g.cm -3 . Comparison between the 
full-statistical (SCO) and SCO-RCG calculations. The maximum number of lines potentially detailed 
per transition array is chosen equal to 800000. 



Figure 7: (Color online) Iron opacity at T=38 eV and p= 4 10 -3 g.cm -3 . Comparison between the 
full-statistical (SCO) and SCO-RCG calculations. The maximum number of lines potentially detailed 
per transition array is chosen equal to 800000. 
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Figure 8: (Color online) Comparison between two SCO-RCG calculations relying respectively on DLA 
and PRTA treatments of lines for transition arrays 3^3/2 —> 5s in a Hg plasma at T=600 eV and p=0.01 
g/cm 3 . The DLA calculation contains 102 675 lines and the PRTA 26 903 lines. 



Figure 9: (Color online) The three independent contributions to photo-excitation calculated by SCO- 
RCG code for an iron plasma at T=193 eV and p=0.58 g.cm -3 (boundary of the convective zone of the 
Sun). 
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Figure 10: (Color online) Number of detailed transition arrays (DLA), partially resolved transition arrays 
(PRTA) and unresolved transition arrays (UTA) for different values of the maximum number of detailed 
lines imposed : 10 2 , 10 3 , 10 4 and 10 6 . For each case, two histograms are displayed: in the first one, the 
detailed calculations are only pure DLA and in the second one they can be either DLA or PRTA. 


a J of configuration C reads 
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(9) 


with ensures that XXjeC = Tc, where Vc is the probability of the genuine configuration given in 
Eq. (0). 

Figure [9] represents the different contributions to opacity (DLA, statistical and PRTA) for an iron 
plasma in conditions corresponding to the boundary of the convective zone of the Sun. We can see that 
the PRTA contribution is of the same order of magnitude here as the statistical one. The calculation was 
performed with a maximum imposed of 10 000 detailed lines per transition arrays. We can see in Fig. |TU] 
that for each value of the maximum number of lines that can be detailed (-/V max ), some UTA are replaced 
by PRTA transition arrays. Of course, the number of remaining UTA decreases with N max . 


4 Interpretation of experimental spectra 

The SCO-RCG code has been successfully compared to several absorption and emission experimental 
spectra, measured in experiments at several laser (Fig. ITTl) or Z-pinch facilities (Fig. U2l) . The comparisons 
show the relevance of the hybrid model and the necessity to carry out detailed calculations instead of full 
statistical calculations. As mentioned in Sec. 0, the quantity which is measured experimentally is the 
transmission, related to the opacity by Beer-Lambert-Bouguer’s law: 

T{hv) = e~P L ^ hv \ (10) 

where L is the thickness of the sample. The relation m between transmission and opacity is valid under 
the assumption that the material is optically thin and that re-absorption processes are neglected. 
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Figure 11: (Color online) Interpretation with SCO-RCG code of the copper spectrum (2 p -+ 3d transi¬ 
tions) measured by Loisel et al. (Loisel et al., 2009, Blenski et al., 2011a & 2011b). The temperature is 
T=16 eV and the density p =5 10 -3 g.cm 3 . 


In SCO-RCG, configuration interaction is limited to electrostatic one between relativistic sub-configurations 
(nij orbitals) belonging to a non-relativistic configuration (ni orbitals), namely “relativistic configura¬ 
tion interaction”. That effect has a strong impact on the ratio of the two relativistic substructures of the 
2 p —> 3d transition on Fig. [TT] 
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Figure 12: (Color online) Interpretation with SCO-RCG code of the iron spectrum (2 p —> 3 d transitions) 
measured by Bailey et al. (2009). The temperature is T=150 eV and the density p=0.058 g.cm 3 . 
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Statistical modeling of Zeeman effect 

5.1 Determination of the moments 

Quantifying the impact of a magnetic field on spectral line shapes is important in astrophysics, in inertial 
confinement fusion (ICF) or for Z-pinclr experiments. Because the line computation becomes even more 
tedious in that case, we propose, in order to avoid the diagonalization of the Zeeman Hamiltonian, to 
describe Zeeman patterns in a statistical way. This is also justified by the fact that in a hot plasma, the 
number of lines is huge, and therefore the number of Zeeman transitions, arising from the splitting of 
spectral lines, is even greater, which makes the coalescence of the spectral features more important. Due 
to the other physical broadening mechanisms, the Zeeman components can not be resolved (Doron et al., 
2014). 

In the presence of a magnetic field B , a level aJ\ (energy E\) splits into 2Ji + 1 states M\ (— Ji < 
Mi < Ji) of energy Ei + gsgiMi , g B being the Bohr magneton and g\ the Lande factor in intermediate 
coupling (provided by RCG routine). Each line splits in three components associated to selection rule 
AM=g, where q =0 for a tt component and ±1 for a a± component. The intensity of a component can 
be characterized by the strength-weighted moments of the energy distribution. The n th — order moment 
reads 


Mn[?]= 3 £ ( _ J M * M ) ( E 2-Ei+» B B{g 2 M 2 - gi Mi}r , (11) 

Mi,M 2 ' 1 q 2 ' 

which can be evaluated analytically (Pain & Gilleron, 2012a & 2012b), using graphical representation of 
Racah algebra or Bernoulli polynomials (Mathys & Stenflo, 1987). 



Figure 13: Effect of a 1 MG magnetic field on triplet transition ls2s 3 S —> ls2 p 3 P of carbon ion C 4+ 
with a convolution width (FWHM) of 0.005 eV. The observation angle 9 is such that cos 2 (0) = 1/3. 
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Figure 14: (Color online) Convolution of Gram-Charlier expansion series with a Lorentzian. The pa¬ 
rameters are a=0.1, <j= 1, 0 : 3=1 and 04 = 5 . The present approach (red curve) relying on a cubic-spline 
representation of the Gaussian and the direct numerical convolution (dashed black curve) are superim¬ 
posed. 
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Table 1: Values of a 3 and ay of the Zeeman components. J< = min(Ji, J 2 ) and J> = max(Ji, J 2 ). 
sgn [x] is the sign of x. 


ak 

e -(k+lf/2 fc 2 (fc + 2)2 + e fc+l /2 ( fc 2 _ X )2 

bk 

e -(fc+i Y/ 2 k{k + 1) [8 + 3k + e fc + 4 / 2 (3 k — 5)] 
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—g -( k + 1 V/ 2 3 + e k+1 / 2 (k — 2) + k 


Table 2: Expression of the coefficients ak, bk, Ck and dk involved in the cubic-spline representation of the 
Gaussian (Eq. (12Tj> ) for k < u < k + 1. 


5.2 Gram-Charlier distribution 

Gram-Charlier expansions are useful to model densities which are deviations from the normal one. The 
expansion is named after the Danish mathematician Jorgen P. Gram (1850-1916) and the Swedish as¬ 
tronomer Carl V. L. Charlier (1862-1934). Historical accounts of the origin of the Gram-Charlier ex¬ 
pansion are given in Hald (2000) and Davies (2005). This expansion, that finds applications in many 
areas including finance (Jondeau & Rockinger, 2001), analytical chemistry (Di Marco & Bombi, 2001), 
spectroscopy (O’Brien, 1992) and astrophysics and cosmology (Blinnikov, 1998) reads 


GC{u) 



_fc=0 



( 12 ) 


where u = (hv — Hi) / y/v, v = ^2 — (Hi) 2 being the variance. The polynomials Hefc(x) can be expressed 
as 


He fc(z) - 2fe /2 H * (^ 2 ) ’ 

where H^.(x) are the usual Hermite polynomials obeying the recurrence relation (Szego, 1939): 


(13) 


(14) 


H fc+ i(x) = 2xH fc (x) - 2fcH fc _i(x) 

initialized with Ho(x)=l and Hi(x) = 2x. The coefficients Ck are given by 

W 2 ] 

j\(k — 2j)\2i 

where [.] denotes the integer part and ak is the dimensionless centered k —order moment of the distribution 
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Table 3: Coefficients 7 Pi i involved in Eq. (l22l) . 


A good representation of the Zeeman profile is obtained using, for each component, the fourth-order 
Gram-Charlier expansion series: 


'i’z(u) 



(o?4 - 3) 
24 


(3 — 6u 2 + M 4 ) , 


(17) 


where 03 (skewness) and ol\ (kurtosis) quantify respectively the asymmetry and the sharpness of the 
component (see Table [T]) (Kendall & Stuart, 1969). This approximate method was shown to provide 
quite a good description (see Fig. [131) °f the effect of a strong magnetic field on spectral lines (Pain & 
Gilleron, 2012a & 2012b). The contribution of a magnetic field to an UTA can be taken into account 
roughly by adding a contribution 2/3 (hbB) 2 ~ 3.35 10 -5 (B(MG )] 2 eV 2 to the statistical variance. 

When all the other broadening mechanisms (statistical, Doppler, ionic Stark) are described by a Gaus¬ 
sian, the resulting profile (convolution of a Gaussian by Gram-Charlier) remains a Gram-Charlier func¬ 
tion with modified moments. However, electron collisional broadening is usually modeled by a Lorentzian 
function 


L(hv, a) 


a 1 
7T a 2 + (hv) 2 


(18) 


as well as natural width. The convolution of a Gaussian by a Lorentzian leads to a Voigt profile (Voigt, 
1912, Matveev, 1972 & 1981, Ida et ah, 2000) but in the presence of a magnetic field, the problem is 
more complicated, since the numerical cost of the direct numerical convolution of a Lorentzian with 
Gram-Charlier function is prohibitive, due to the huge number of lines involved in the computation. It 
reads 


/ OO 

GC (u) L(t - u, A )du 

-OO 

where t = hv/a. a = yjv being the standard deviation of the distribution and A = a/a. 


(19) 


5.3 Convolution of a Lorentzian function with Gram-Charlier expansion se¬ 
ries 

The convolution product (1191) requires the evaluation of a cumbersome integral which reads 


C(t) 
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In order to fasten the calculation, the Gaussian is sampled at the points u = —to, — m+ 1, • • ■ , 0, • • • , to — 
1, to (in practice we use m=6 ) and interpolated using cubic splines (de Boor, 1978) on each interval 
[ k , k + 1] by the formula 


e “ 2 / 2 = afc + bk u + Ck u 2 + dk u 3 . (21) 

The coefficients a,k, bk, Ck and dk in the interval [k, k + 1] are determined by the continuity of the function 
and its derivative at the points u = k and u = k+ 1. The resulting expressions are given in Table [2] The 
Gaussian is assumed to be zero for |u| > to. Limiting Gram-Charlier expansion series to fourth order 
(see Eq. (1171) 1. one now has to deal with the convolution of a Lorentzian by a polynomial of order 7. This 
can be writtero 

m— 1 7 

t p.fc Sp,k(t) , ( 22 ) 

p=—m fc=0 

where the coefficients 7 Pi fc of the polynomial are given in Table [31 and 


7TV 27T( 


r 

Sp, k (t) = / 

J V 


P + 1 


A 2 + (t - m ) 2 


r C?M 


= )e(‘) [* (“w) - (¥)] ■ 

^=0 A / 


(23) 


The function ( pe(w ) is equal to 
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,«+i 
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w ) = 2 E 1 ( 4 3 2 ;-«>), 


(24) 


where 2 E 1 is a hypergeometric function, but can be efficiently obtained using the recurrence relation 


<j>e(w) = u> e 2 - (t>e- 2 (w) (25) 

with 

<j> o(w) = arctan(w) and (j)\{w) = — In [l + w 2 ] . (26) 

Such a method provides fast and accurate results, even for very asymmetrical and sharp Gram-Charlier 
distribution (see Fig. fl4l) . The total line profile results from the convolution of ^z with other broadening 
mechanisms. If cr < a/10, we take only the Lorentzian L(hv,ci). On the other hand, if a < er/150, we 
keep the Gaussian. If one has to convolve C{t) by an additional Gaussian of variance a' (representing 
Doppler broadening for instance), cr, a 3 and oq must be replaced respectively by: 

( G — V <7 2 + <J /2 

l b 3 =a 3 (|)' (27) 

l «4 = a 4 (f) 4 

Figs. fl5l and UHl illustrate the impact of a 10 MG magnetic field (typical of ICF) in the XUV range for a 
carbon plasma at T=50 eV and p=0.01 g/cm 3 . 

2 In the general case, the upper bound of the sum over k is equal to JV+3, where N is the order of the Gram-Charlier 

expansion series. 
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Figure 15: (Color online) SCO-RCG calculations (transitions Is —► 2 p) with and without magnetic field 
for a carbon plasma at T =50 eV and p=10 -2 g.cm -3 (conditions typical to ICF). 
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Figure 16: (Color online) SCO-RCG calculations (transitions Is —> 2 p) with and without magnetic field 
for a carbon plasma at T= 50 eV and p=10 -2 g.cm -3 (condidions typical to ICF) in a spectral range 
close to the one of Fig. [15] 
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6 Conclusion 


By combining different degrees of approximation of the atomic structure (levels, configurations and su¬ 
perconfigurations), the SCO-RCG code allows us to explore a wide range of applications, such as the 
calculation of Rosseland means, the generation of opacity tables, or the spectroscopic interpretation of 
high-resolution spectra. The Partially Resolved Transition Array model was recently adapted to the hy¬ 
brid statistical / detailed approach in order to reduce the statistical part and speed up the calculations. 
An approximate approach providing a fast and quite accurate estimate of the effect of an intense magnetic 
held on opacity was also implemented. The formalism requires the moments of the Zeeman components 
of a line aJ —> o! J', which can be obtained analytically in terms of the quantum numbers and Lande 
factors and the profile is modeled by the fourth-order A-type Gram-Charlier expansion series. We also 
proposed a fast and accurate method to perform the convolution of this Gram-Charlier series with a 
Lorentzian function. Such an algorithm is useful in order to account for distorsions of the Voigt profile, 
since the direct numerical evaluation of the integral becomes rapidly prohibitive. More generally, It can 
be helpful for models relying on the theory of moments (Bancewicz & Karwowski, 1987), used in most 
opacity and emissivity codes. In the future, we plan to extend the statistical modeling of Zeeman effect 
using temperature-dependent moments (see Appendix) and to improve the treatment of Stark broadening 
in order to increase the capability of the code as concerns K-shell spectroscopy. 
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A Temperature-dependent moments of Zeeman Hamiltonian 

Polarized synchrotron radiation can be used to determine the magnitude, the orientation and the tem¬ 
perature and magnetic-field dependence of the local rare-earth magnetic moment in magnetically ordered 
materials. Thole et al. (1985) proposed a theory which predicts an anomalously large magnetic dichroi'sm 
in the X-ray absorption-edge structure. The square of the matrix element of an optical dipole transi¬ 
tion from a state aJM to a final state a'J'M' is, according to the Wigner-Eckart theorem, proportional to 
the square of the 3j symbol times the reduced matrix element (line strength in the absence of a magnetic 
field): 


S aJM , a ’j'M> = ( \ Jm> ) \(<* J \\ cil) \W J ')\\: ( 28 ) 

where a labels different levels of equal J and g=0 for light polarized in the field direction, q = ±1 for right- 
or left-circularly polarized light perpendicular to the field direction. The partition function associated to 
a particular level aJ reads 


e -C aJ M _ sinli [ C a j(J + 1/2)] 
“ , ’ sinh(C a j/2) 


(29) 


where C a j = gBBg a j/(ksT), fig being the Bohr magneton, B the intensity of the magnetic field and 
g a j the Lande factor of level aJ in intermediate coupling. The first order moment can be expressed as 


( M ) = J- ^2 Me~ CajM = —(J + 1/2) coth [C aJ (J + 1/2)] + i coth(C , Q ,j/2) 

aJ M— — J 

= —J Bj(C aJ J), (30) 

where Bj denotes Brillouin’s function (Darby, 1967, Subranranian, 1986): 

.. 2 J +1 /2 J +1 \ 1 f x \ , . 

BjW = ^^ co,h (^^7“w coth ( w)- (31) 

For (M 2 ), one has 

1 J 

(M 2 ) = — M 2 e~ c “ jM = J{J + 1) + (M) coth(C a j/2) (32) 

aJ M=—J 

and the higlrer-order moments can be obtained using the following relation: 

/ -| \ yl CSTL ry O 

(M n ) = dc<r ^ = {M)(M n ~ 1 ) - —(M- 1 ), (33) 

where n £ N. 
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